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Among the mechanisms for lattice structural deformation, the electron-phonon interaction me- 
diated Peierls charge-density-wave (CDW) instability in single band low-dimensional systems is 
perhaps the most ubiquitous. The standard mean-field picture predicts that the CDW transition 
occurs at all fillings and all values of the electron-phonon coupling g and the adiabaticity parameter 
t/u!o- Here, we correct the mean-field expression for the Peierls instability condition by showing 
that the non-interacting static susceptibility, at twice the Fermi momentum, should be replaced by 
the dynamic one. We derive the Luttinger liquid (LL) to CDW transition condition, exact to second 
order in a novel blocked perturbative approach, for the spinless one-dimensional Holstein model in 
the adiabatic regime. The small parameter is the ratio gujo/t. We present the phase diagram at 
non-half-filling by obtaining the surprising result that the CDW occurs in a more restrictive region 
of a two parameter {g'^ujo/t and t/ujo) space than at half-filling. 

PACS numbers: 71.38.-k, 71.45.Lr, 71.30.-hh, 75.10.-b 



I. INTRODUCTION 



Over the last few decades, electron-phonon interac- 
tion physics has offered a variety of intriguing and ex- 
citing phenomena such as superconductivity (inorganic 
and organic), CDW states, colossal magnetoresistance, 
metal-insulator transition, polaronic ordered phases, etc 
[1]. Of the electron-phonon models available, the spin- 
less Holstein model [2,3] is a simple and widely used 
model that mimics strongly correlated electron systems 
with a strong on-site coulombic repulsion and short range 
electron-phonon interactions. It was conjectured a long 
while ago by Peierls that ID electron-phonon metallic 
systems, such as the spinless Holstein model, will undergo 
an electronic charge-density-wave (CDW) transition with 
a concomitant lattice distortion of the same periodic- 
ity [4]. Quasi- ID organic charge transfer salts [such as 
TTF(TCNQ)] and conjugated polymers [such as (CH)^] 
as well as inorganic blue bronzes (e.g., K0.3M0O3), mixed 
valence Platinum chain compounds (e.g., Krogmann's 
salt), and transition metal chalcogenides (e.g., NbScs) 
[5-8] exhibit such spontaneous symmetry breaking in the 
ground state and are good candidates for the Holstein 
model. Furthermore, even the strongly correlated two- 
band manganite systems (e.g., Lai-xCaxMnOs) [9], in 
the low-doped regime, can be modeled using the Holstein 
model [10,11]. 

Contrary to the mean-field picture, at half-filling, it is 
now clear that the LL to CDW transition occurs only 
above a critical electron-phonon coupling strength that 
depends on the adiabaticity. Progress has been made 
over the last few decades, in terms of studying the Peierls 
transition at half-filling of the spinless Holstein model, 
by using various techniques such as quantum Monte 
Carlo simulations [12-15], two-cutoff renormalization- 



group analysis [16], variational method [17,18], density- 
matrix renormalization group (DMRG) method [19], and 
exact diagonalization [20]. However, a controlled ana- 
lytic treatment of the quantum phase transition has been 
reported only recently and that too only in the anti- 
adiabatic regime [21]. Contrastingly only little effort has 
been devoted to understand the quantum phase transi- 
tion away from half-filling [22]. 

The present paper is aimed at providing a well con- 
trolled analytic approach to understand the Peierls quan- 
tum phase transition in the adiabatic regime for the one- 
dimensional spinless Holstein model at a general filling. 
We employ a novel blocking approach that avoids the dif- 
ficulties posed by both time-dependent- and degenerate- 
perturbation theories. Using the condition that the ef- 
fective phonon frequency becomes soft at the symme- 
try breaking point, we obtain an instability criterion 
that preempts the energy levels crossing condition for 
phase transition. In the adiabatic regime and at half- 
filling, we capture the essential features of the LL-CDW 
transition results obtained by the "benchmark" DMRG 
method in Ref. [19]. At fillings other than half, we show 
that the LL phase certainly exists in the small polaron 
limit when g'^ujo/t >> max(l, i/wo) and in the extreme 
anti-adiabatic regime. Furthermore, in the adiabatic 
regime, we also demonstrate that CDW phase does exist 
at intermediate values of the electron-phonon couplings 
g < t/ujQ. We propose a qualitative phase diagram as a 
guide for future work. 



II. PHONON SOFTENING IN THE LL PHASE 

We begin by considering the phonon-softening in the 
LL phase as the signal for lattice deformation in the ID 
Holstein model. The non-interacting Hamiltonian 
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and the perturbation 



(1) 



(2) 



together make up the Holstein Hamiltonian. In the above 
equations, = ct,_^^cg, is the density operator, eg, 

is the electron destruction operator with k' Umited to 
the first Brillouin zone, e^, = — 2icos(fc') with lattice 
constant being taken to be unity, t is the hopping in- 
tegral, Cq- is the phonon destruction operator, luq is the 
optical phonon frequency, and N is the number of sites. 
The eigen states and eigen energies of Hq are given by 
= \n;m) = \n)ei\m)ph (with l^o) = |0;0) being the 
ground state with zero phonons) and E^^ respectively. 
Whereas, for the interacting Hamiltonian H = Hq + Hi, 
the corresponding eigen states are and the eigen en- 
ergies are . 

Now the double time derivative of an operator A is 
given by 



A^~[[A,H],H]. 



(3) 



From the above equation, when A is taken to be the 
ionic position coordinate 



^l/{2Mu;o){ap + al^^) 



and upon making the static mean-field approximation 
p-fj oc Xo{p, 0)Qp [8], we obtain the following expression: 

Qp = -u^ [1 + 2g'^u;oXo{p, 0)]Qp = -uJrQp- (4) 

We know that in ID the non-interacting polarizability 

Xo{p, 0) has a negative divergence at wavevectorp = 2kF- 
Thus it appears that the renormalized phonon frequency 
ojR becomes soft even for vanishingly small electron- 
phonon interaction leading to lattice distortion. We will 
now proceed to derive the true phonon softening con- 
dition rigorously. To this end we calculate the matrix 
elements of Eq. (3) and obtain 



(5) 



When = {E^^ — E^^Y < 0, instability occurs for 

transition from 1$,,) to |$„i) provided that ($„, |v4|$„) ^ 
0. For the total Hamiltonian Hq + XHi, we obtain 
perturbatively to be E^^ = E^^ + X'^E^^ + ... where 

E^^ is the second order correction to the energy. Then, 
to second order in A, we obtain 



= (^L - El? + 1X\El^ - ElME^l - (6) 



Thus to order in perturbation, as the strength of the 
interaction is increased, = before _E<i,,,^ = . 

We will now focus on the eigen states and the eigen 
energies of Hq + Hi to obtain the phonon softening 



condition. The energy of the (expected) ground state, 
upto second order in perturbation, is given by: 



-£'*n — Tq 



9 ^0 
N 



\{m\elPif\0)el\ 



(7) 



where Tq ( = E^^) is the non-interacting kinetic energy 
of the l^o) state and = ?m — Co with being the 
energy of |m)e(. Let \'ip'^) = |0;n_p) with \n_p)ph cor- 
responding to a state with n phonons all of which being 
in the —p state. Then, the corresponding interacting 
state 1 5"") yields the energy difference Eq,o — E^s,o^ = 
wq + E(— p, Wo) where the self-energy, which is complex 
in general, is given by T,{p,ujq) = g^u!^Xo{p,(^o) [see Ap- 
pendix A for details] with the non-interacting (Lindhard) 
polarizability xo being defined as [23] 



Xo{p,uJo) = -^Yl 



{m\elpp\0)elf \{m\elP-p\0)el\ 



Wo - CmO + Wo + CmO + iV 



Hence we see that, although the above energy difference 
yields the expression = loq[1 + 2g'^uioXo{Pj'^o)] [based 
on Eqs. (5) and (6)] with the form being similar to that 
in Eq. (4), the complex nature of the self-energy com- 
plicates identifying the phonon softening condition for 
lattice instability [24]. We will adopt an alternate per- 
turbative procedure to obtain the lattice instability cri- 
terion. 



III. HALF-FILLED CASE 

We begin by observing that the expression 

Wo = ^k-p ~ ^fe = 4i sin(A;ir - fc) sin(fcir), (8) 

with p = 2kF has two solutions for k. Then at half- 
filling [kp = 7r/2) and k < kp, the two solutions to 
Eq. (8) are ±fc = tt/2 — arcsin[a)o/4t]. Thus, the 
states 1^}) = jfc — » k — p;0) (corresponding to excit- 
ing, from the ground state, the electron at fc to fc — 
and IV'i) = I — fc ^ —k — p;0) are degenerate with 
\tpi) = |0, l-p-) and are connected to \tpi) through Hi 
as {tpi l-ffilV'i) 7^ 0. Hence, employing degenerate per- 
turbation theory seems to be a natural choice to study 
lattice period doubling. However, the number of degen- 
erate states that need to be considered increases linearly 
with number of phonons in the state |?/;") = |0; n_p) [see 
Appendix B.l]. Then, to calculate for large n becomes 
difficult! 



2 



.-X--. 



✓ N 
/ \ 



J \ L 



k-p -kp -To k 



-k-p 



occurring ai g > gs pre-empted by the phonon soften- 
ing occurring at g > gc = gs/V^ with gc defined by the 
following expression: 



1 + 2glujQRexo{p,(^o) = 0. 



(12) 



FIG. 1. Blocking procedure at half-filling. Electrons 
at Tkp ~ T'i"/2 are excited from the ground state (dark line) 
to ±fc — p to block the excitations ±fc — > ±fc — p respectively. 

To circumvent the above problem, we adopt the fol- 
lowing approach. Wc choose our starting state to be 
\(j)n) = \ — kp ^ k — p, kp ^ —k — p; n-^) (n > 1) such 
that, by exciting the electrons at T^f from the ground 
state to ±fc — p outside the Fermi sea (FS), the excita- 
tions k ^ k — p and —k^ —k — p are Pauli blocked [as 
shown in Fig. 1]. Now, ikp have been excited for ease of 
mathematical manipulation. The state upon turning 
on interactions, yields the energy [see Appendix B.2 for 
details] : 

Eg>„ = E^„ + in + l)wo[l -I- g'^LOaRexoip, loq)] 

+g'^(jJoRe[2Uo{k - p,kF,(^o) - Xo{p,oJo)], (9) 

where 

ReHo (s, Q,^o) = --^^ 



1 — 



Tip 



ep + LJo — eg ojo — ef 



es + ojQ — ef ef + u)o — ^3 



For Wo/(4f) « 1, 
Re[2no(fc -p, fcF,wo) 



Xo(p,Wo)] ~ 



8t 



(10) 



The above approximation underestimates the actual 
value of 21X0 — Xo by less than 5% for t/wo > 4. In Ap- 
pendix C, Fig. 6(c) shows that 2no — xo > for f/cjo > 1. 
In the above Eq. (9), for n 00, energy instability 
(i5$„ — E^g < 0) occurs for values of g larger than gs 
given by 1 + .g|cjoRcxo(p, ^o) = 0. For g > gs, En has 
no lower bound which is an unphysical situation. Also 
when g > gs, -B*„+i - -B$„ = 1-1- .g^woRexo(p, ^^o) < 
for all n > 1 which leads to the remarkable situation that 
all cross at the same g = gs- To second order in the 
small parameter of perturbation, similar to Eq.(6), one 
obtains 



{E^„-E^„f « (n + l) 



+2(n+l)^ln 



25^woRexo(p,wo)] 



9ujo 



(11) 



Thus we see from Eq. (11) that, in the adiabatic regime 
and for large n, the above mentioned energy instability 



The above equation is one of our main results and is 
the correction to the mean-field instability condition 1 -|- 
2g^u}oXo{'P,0) = obtained from Eq. (4). The operator 
A, which produces non-vanishing matrix elements in Eq. 
(5), is given by A = 4_/-^/2cl^fe_/^/2(aLp-)"- Thus, 
the system becomes unstable towards absorbing a large 
number (n) of phonons leading to a macroscopic deforma- 
tion as explained below. The displacement-displacement 
correlation function is given by 



{'S>n\QlQj\^n) = (KlQlQj 



_ ncos[{j-l)p] 



Thus we see that, for non- vanishing values of n/N, one 
obtains an observable ionic-position modulation when 
is the interacting ground state. The above Eq. (13) 
is true for all fillings and for any eigen state with n_p- 
phonons. 
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FIG. 2. Critical coupling gc versus adiabaticity pa- 
rameter t/wQ. Comparison of the gc values obtained at (a) 
half-filling in this work, in Ref. [19] using DMRG, and in Ref. 
[16] using two-cutoff RG; and (b) various fillings (y) in this 
work. 



3 



The critical coupling cjc. given by Eq. (12), can be 
expressed analytically as follows: 



TT 



In 



(1_ 0:3^)2 tan 

(l + 0^^)2-(7tanfcF)2 



-2 In 



1 + 



(14) 



where 7 = a;o/(4t sin(A;F)) < 1- Fig. 2 depicts, for vari- 
ous filling factors the variation of the critical coupling 
Qc with the adiabaticity parameter t/uja- At half- filling, 
for values of t/wo > 5, our theoretical curve is quite close 
to the numerically determined values of gc as reported in 
Ref. [19]. Furthermore, at half-filling and for 7^ << 1, 
our expression for g^. [given by Eq. (14)] reduces to the 
two-cutoff renormalization result of Caron and Bourbon- 
nais [16], i.e., = 2rfexp(— Trt/g^wo) (with c ^ 1), when 
we take c = 4. The numerical agreement between the two 
expressions is depicted in Fig. 2(a). 



IV. LESS THAN HALF-FILLING 



We will now consider fillings that are less than half- 
filhng. The line depicted by Wo = l/[4sin(2fc_F) sin(fc_F)] 
[obtained by setting k = —kp in Eq. (8)], corresponds 
to the divergence of Xo{P^^o)- In the region above 
t/wo = l/[4sin(2A;ir) sin(A;i;')] (see Fig. 3), the excita- 
tion energy expression - — eg = with p = 2kF is 
satisfied by one wavevector for < kp [as seen from Eq. 
(8)] (see Appendix C for a complete analysis). We will 
now consider the region t/wo > l/[4 sin(2fc_F) sin(A;F)]- In 
this region, the degenerate states are [?/'") = |0; and 

Itpli) = |fc — > k — p;{n — We obtain the lower 

eigen energy, which corresponds to the state |'^~) = 
[IV'") - IV'i)]/\/2, to be (see Appendix C) 

E^- = +E'^ + nwo[l + .g^a)oRc=Xo(p, wo)] 

-0.5/w^Re[xo(p,wo) - no(fc - j3,fc,wo)], (15) 




FIG. 3. Curves relevant for identifying different per- 
turbative regimes. Plot of t/cjo and filling factor i' = kp/Tr 
values satisfying Rexo(2fcF,wo) = and curves for the func- 
tions t/ioo ~ l/[4sin(A;F) sin(2A;F)], = l/[4sin(A;F)], and 
= l/[2sin(2fcF)]. 

It should be pointed out that only the region 
Rexo(l5, t^o) < is relevant in obtaining gc- The curve 
Rexo(p, t<^o) = is depicted in Fig. 3 and exists only 
for kp < 7r/4. Above (below) this curve, Rexo(Pj'^o) is 
always negative (positive). Furthermore, for kp > 7r/4 
and all values of the adiabaticity parameter t/uJo above 
(below) the hue t/tJo = l/[4sin(fei^)], one can show ana- 
lytically that Rexo(p, t^o) is always negative (positive). 

It can be shown that, for a given filling v and any value 
oit/ojQ where Rexo(P!"^o) < 0, the macroscopic instabil- 
ity condition is always given by Eq. (12) (see Appendix 
C). For filling factors above 0.25, as shown in Fig. 2, 
the Qc decreases with decreasing t/uo with a downward 
kink appearing at a certain value of t/uo corresponding 
to the negative divergence of xo(P;"^o)- At the point of 
divergence of xo{p,^o), perturbation theory is no longer 
valid. For kp < n/A, the gc initially decreases with de- 
creasing t/uJo until a certain value of t/uio', while below 
this value t/ojo, the value of gc again increases due to the 
fact that Rexo(p, 1^0) value approaches zero value (see 
Fig. 3). Lastly, we would like to point out that the gc 
values are not reliable when gc(^o /t > I and hence in the 
entire anti-adiabatic regime {t/ojo < 1) the gc values are 
suspect [see Fig. 2(b)]. 



where = —^=goJo is the first order energy correc- 
tion. In arriving at the above energy, we have ignored the 
contribution 'T-/[A^(e^_2p ~ "-^o)] corresponding to 

exciting the electron at k—p to the state k—2p by destroy- 
ing a phonon of momentum —p. This is valid provided 
+ 0. The case when eg_2--eg_, -Wo = 
(i.e., i/wo = l/[2sin(2A:j?)]), will be discussed in the Ap- 
pendix C.l. Obviously, the macroscopic deformation in- 
stability condition is still g > gc with gc given by Eq. 
(12). 



V. T=0 PHASE DIAGRAM AT 
NON-HALF-FILLING 

At non-half filling, we will now discuss the quantum 
phase transition based on the perturbation theory re- 
sults derived above and the work reported in Ref. [21]. 
In the extreme small polaron regime, for (7 > 1, it was 
shown earlier that the effective Holstein Hamiltonian can 
be recast as an effective spin Hamiltonian (using Wigner- 
Jordan transformation) as follows (see [21] for details): 
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j 3 



, (16) 



where C, = t/g^ujQ is the polaron size parameter. 
The above; equation was obtained by assnining that 
the phonons are frozen in the Lang-Firsov trans- 
formed (LFT) phononic ground state exp[5Ej c]cj(aj — 
aj+i)]|0)ph [25]. Now, for the above Eq. (16) to be the 
basis for studying phase transition, each of the coeffi- 
cients of the second, third, and fourth terms on the rhs 
should be significantly smaller than ljq (so that the LFT 
phononic ground state remains unaffected). In Eq. (16), 
for Q « 1, the coefficients of the nearest neighbor and 
the next to nearest neighbor interactions in the trans- 
verse direction are much smaller than wq- Contrastingly, 
the coefficient of the nearest neighbor interaction in the 
longitudinal direction is much smaller than lhq always 
when t/ujo < 1; while for t/u)o > 1, it is much smaller 
only when [t/{gu)o)]'^ « 1 (i.e., for large values of g). 
Note that, when C << 1, the last term is negligible. Us- 
ing Bethe ansatz, we know that anisotropic Heisenberg 
model always yields a Luttinger liquid away from half- 
filling [26-28]! From the above analysis, it follows that 
a LL results for all values of l/( » max(l, f/wo). Fur- 
thermore, for t/uo « 1 and g > 1, the above Eq. (16) 
is valid with the last term on the rhs being negligible and 
consequently LL results away from half-filling. 

Next, when g < 1 and t/uo « 1, we get the cor- 
responding effective spin Hamiltonian from the effective 
Holstein Hamiltonian to be [21] 



E^i+ce-^'E(4^7+i+H-c.) 

j j 



(17) 



In arriving at the above equation too, it was assumed 
that the phonons are in the LFT phononic ground state. 
Such an assumption is justified because the coefficients 
of the second, third, and fourth terms on the rhs of the 
above equation, are much smaller than wq- In Eq. (17), 
the small parameter is t/ujo and the last two terms are 
negligible compared to the second term when the adia- 
baticity parameter t/uio « 1- Then, this implies a LL 



state for all fillings. For and any value of t/wo, we 

do not expect a CDW state. 

The phase diagram (see Fig. 4) is drawn qualitatively 
for a general filling away from half-filling. For v < 0.25, 
the CDW region shifts to the right with decreasing i/ as 
can be surmised from the region of validity gc^^o /t < 1 in 
Fig. 2(b). The regions where LL is certain is indicated. 
For t/ojo > 1, since we need [t/{givo)]^ « 1 for the va- 
lidity of Eq. (16), the boundary of the LL-certain-region 
is linear and of the form g^u^jt = Dt/ujQ where the slope 
D » 1. Furthermore, for guio/t << 1, one expects a LL 
phase and hence we get a linear boundary (of the form 
g'^ujo/t = dt/ujo with d << 1) for the LL phase in the 
lower left part of the diagram. 
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FIG. 4. Zero temperature phase diagram at non-half 
filling. Regions where CDW and LL phases certainly exist 
are depicted. The calculated transition from LL to CDW is 
indicated by a dark line along with its expression. 

The thatched portion corresponds to a region where 
a CDW state is certain. The upper boundary of the 
thatched region, which is obtained by the condition 
givo/t ~ 1 and is therefore linear with slope of order 
unity, corresponds to the breakdown of the perturbation 
theory used in identifying the CDW transition. The only 
certain location of the transition from LL to CDW is in- 
dicated by a solid dark line and is approximately given 
by g'^LUo/t ^ l/(ln(8i/a;o) [while the exact relation is ex- 
pressed in Eq. (14)]. However, it is unclear in the rest 
of the figure where exactly the transition from a CDW 
state to a LL state takes place. 

Finally, it should be emphasized that there are two 
different types of phase transitions. The nature of the 
phase transition on the adiabatic side, at intermediate 
values of g < t/ojo, is driven by a macroscopic ionic lat- 
tice distortion. Quite differently, in the restricted region 
of the small polaron limit where 1/C >> max(l, t/cjo) 
and in the extreme anti-adiabatic regime {t/wo « 1), 
the CDW is driven by a small-polaron-interaction based 
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mechanism. In the latter case, the coordinate Bethe 
ansatz argument precludes the possibility of a small po- 
laronic CDW away from half-filling in sharp contrast to 
the half-filled case [19,21]. 
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APPENDIX A 



The interacting state 1^1/°), corresponding to the non-interacting state |^°) = |0;n_p-), yields the following energy 
expression: 



£^^0 = To + nujo — ^ 



<t„^^o L <l>i n 



(Al) 



where Tq = -E? = —2Ntsm{kF)/T^ is the non-interacting kinetic energy of the ground state |(/)o) and 



E 



2 2 

^5 Wo 



J_ \{m\elPq\0)el\''' 

N J-'^ Uo + Wo 



»^Xo(p,wo) 



(A2) 



with T] 0+. Then the above Eq. (A2) yields the energy difference E^yo ^ — Eq,a = ljq + E(— p, wo) where the 
self-energy of a phonon S(p, Wq) = Q^^oXoiPti^o)- The self-energy is displayed in Fig. 5 with the bubble representing 
the polarizability Xo{P: ^o) and each of the electron-phonon interaction vertices corresponding to the factor giuo- 



q + p, 8+0) 



P' (0, 




p. (0, 



FIG. 5. Self-energy of a phonon. The solid and the dashed lines depict the electron and the phonon propagators respectively. 



APPENDIX B 



1. The set of degenerate states for half-filled case 



The solutions of the expression 



Wq = ec_ir~ ^fc = 4tsin(fcir - k) sin(fcir), 



(Bl) 



for p = 2kF are given by 



k = kp — arcsin 



Wo 



4f sin(fcF) 



(B2) 



and 



k = —IT + kp + arcsin 



Wo 



At sin(fcF) 



(B3) 
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Then, the set of states that arc degenerate with |0;n_jj) and that should be considered in degenerate perturbation 
theory consists of the foUowing states |0; {n — '2.m) ^.jfm^j^^jfn , \k — > k — p;{n — 2m — ^)~frn^j:_pm_^i^_J), 

I - A; -k - p;{n - 2m - l)-pm^^_pm_^^_p), \k ^ k - p,-k ^ -k - p; {n - 2m - 2)_ffm^^_pm_^j:_p), \k 

-k-p; (n- 2rn- 2)_j;(m+ l)2^_p.m_2£_-), \ - k ^ k - p; {n - 2m - 2) _jfln^j:_Jm + 1) , where m=0,l,2,3,... 

with the constraint that the number of phonons is non-negative. Thus we see that the number of degenerate states 
increases Hnearly with n and is given by 3n. 



2. Derivation of energy in the half-filled case 



The starting state \(pn) 
the following energy: 



kp ^ k — p,kF ^ —k — p; ri-p) (n > 1), after switching on the interactions, results in 

'mHA<t>n)?' 



= To + (n + l)uJo - 



E2. -E° 



(B4) 



To evaluate the last term on the right hand side (rhs) of Eq. (B4), we use Eq. (A2) and obtain 



m\HM\' 



E 



Re 



E° -EO 



where 



Reno(s, 9> t^o) = --^ X! 



1 — np 



g a;oRe[no(fc - p, -kp, u>o) 



+ no(-fc -p,kF,ujo)] 
-2g^ujQ Rello {k - p,kF,(^o), 



(B5) 



ep + loq — eg eg- + Wo — Cf 



+ t^o — Er ep + UlQ — eg 



(B6) 



In the above expression for Ho, the first and second terms correspond to adding the contributions due to the electron 
at s going to a state outside the Fermi surface (FS) and those due to the electrons within the FS going to the state 
q respectively. Whereas the third and fourth terms, on the rhs of Eq. (B6), represent subtracting contributions due 
to electrons within the FS going to the state s and those due to electron at q going outside the FS respectively. In 
obtaining Eq. (B5), the terms that are ignored or overcounted are negligible for large N. Prom the main text we 
know that the state \4'q), upon turning on the interaction, yields the energy 



To 



5(2^2 ^ \{mUpq\G)el\^ 



10 \^ 
r 2-^ 



(B7) 



Then, from Eqs. (B4)-(B7), we obtain 

= + (n + l)wo[l + g'^uJoRexoiv^ ^o)] 
+52a;2Re[2no(fc - p, kp, wq) - xo{p, wq)]. 



(B8) 



In the above equation, 2no — Xo > for all values oi t/ojo > 1. 
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APPENDIX C 




FIG. 6. Plots of Qc and g\ (drawn till their crossing point) versus t/ujo when the excitation expression eg - — eg = wo has 
only one solution k (with k < = p/2) for (a) u = 0.1 and (b) v = 0.3 and two solutions for (c) v = 0.5. 

Here, we consider in detail the non- half- filled case. For kp larger (smaller) than 7r/4, the line given hy t/wo = 

l/[4 sin(2fci?) sin(A:F)] corresponds to the larger (smaller) magnitude wavevector solution k lying on the Fermi surface 
as seen from Eq. (B3) (Eq. (B2)). In the region above (below) t/ioo = l/[4 sin(2A;i;') sin(fci?)] (see Fig. 3 in the 
main text), the excitation energy expression = Wo with p = 2kF is satisfied by one (two) wavevector (s) for 

< kp > 7r/4 [as seen from Eq. (B3)] whereas for < kp < 7r/4 it is satisfied by one (zero) wavevector(s) [as seen 
fromEq. (B2)]. 

In the region t/cjo > l/[4sin(2A:i?) sin(fc_F)], the degenerate states are 1-0°) = |0;n_p') and ji/'i) = \k^k—p;{n — 
Unlike at half-filling, only two (i.e., the above given two) degenerate states need be considered for carrying out 
degenerate perturbation theory. The two basis states are \tp:^) = [\i^n) ~ li'n)]/'^^ I'^n) = [IV'n) + IV'n)]/v^ with 
{\tpn}) yielding the lower (higher) eigen energy. Then, from degenerate perturbation theory, one gets 



E.,=E,,-^g.,- ^ . (CI) 



where the second term on the rhs is the first order energy correction; furthermore, it is understood that the states |(/);) 
do not belong to the subspace spanned by jV'n )• Next, to evaluate the last term on the rhs of the above Eq. (CI), we 
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use the following non- mixing fact: 



E 



2 E E° -E° 



Using Eq. (A2) and on noting that 



2 ^ El -E° 



(C2) 



E 



E" 



- E 



\{m\elPg\0)e 



-(n - l)Rexo(p, wo) - Reno(fc - p, fc, wo) , 



we obtain from Eqs. (C1)-(C3) 



E^- = E<s,^ +E'^ + na;o[l + g'^uJo^Xo{p, wo)] 

-0.5fl'^WoRe[xo(p,wo) - no(fc -p,fc,wo)]• 



(C3) 

(C4) 
(C5) 



In Eq. (C3), ILo{k — p, k,LOo) has been obtained in a manner similar to that in the half filled case. As mentioned in 
the main text, the macroscopic instability condition is given by 



1 + 2gluJoRexoiP, ^o) = 0. 



(C6) 



However, for each filling below half-filling, there is a corresponding critical value of the adiabaticity parameter {t/uo) 
below which Re[xo — Hq] > in the above Eq. (C5). Consequently, the interacting state I'i'i) produced by the 
one phonon state \tpi) is more stable than the interacting state |$o) a-t an clcctron-phonon coupling value given by 
gi < g < gc where gi corresponds to {E^- — = [sec Fig. 6]. Only above this critical adiabaticity parameter 

value do we have l^*^) , for n —* oo, as the most stable state at a coupling g > gc < gi- Furthermore, it should also 
be noted that the interacting state is also a LL. Thus, it is obvious that the LL to CDW transition occurs only 
at g = g'^ . Lastly it should also be mentioned that, instead of using the above degenerate perturbation theory to 
obtain the instability condition, one can also adopt a blocking approach similar to that at half-filling by exciting one 
electron at the Fermi surface to the state k — p and then employ non-degenerate perturbation theory. In the latter 
case, the condition for macroscopic instability is still the same while the critical adiabaticity parameter value (above 
which gc < gi) is slightly larger. 



1. The case when e^_p - eg = eg_,^^ - = jjp 

In the non-half-filled regime, we will now consider the special case where the excitation expression - — eg = wq 
and £j^_2j^^ ^'k-p ~ ^0 E^re simultaneously satisfied by one k (with k < kp = p/2). This will hold when 2t cos(A') = ujq 
and consequently when t/ojo = l/[2 sin(2A;F)]. For values of the adiabaticity parameter t/coo and kp that lie on the 
line t/u)o = l/[2sin(2fcir)] depicted in Fig. 3 of main text, the method involving only two degenerate states breaks 
down. To analyze the CDW instability, one can use a blocking method similar to that used at half-filling in the 
main text. We Pauli block the states k — p and A: — 2p by the two electrons on the Fermi surface. The blocked state 
= Ikp k — p, —kp k — 2p; ri-^) yields the energy for the interacting state |$^) to be 

Eis>i^ = E^„ -l-a;o(l + csc(fcF)) + nwo[l + g^woRexo(p, wo)] 

+g^ujlRc[Uo{k -p,kF,LUo)+Uo{k- 2p, -kp^ujo)]. (C7) 

Here too the critical coupling gc, for macroscopic instability, is still given by Eq. (C6). However, for filling factors less 
than approximately 0.22, the one phonon interacting state \^[) is the lowest energy state for gi < g < gc (see Fig. 
7). At fillings above 0.22, the large n phonon interacting state is the most stable state for g > g,, < gi. It should 
also be mentioned that the above blocking procedure can also be used when e^r^ - — = uq and eg_2p' ~ ^k-^ ~ '^o in 
which case the values of t/u)o and kp lie close to the curve t/oJo = l/[2 sin(2A:i?)]. 
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FIG. 7. Plot of Qc and g\ versus fillings v when the excitation expressions eg_^ — ejj = wo and ^%_2^ — = '^o are 
simultaneously satisfied by one k (with k < kp = p/2). 

2. When both the solutions exist 

Here we will consider the lattice instability for the case i/wo < l/[4sin(2A:F) sin(A:F)] when two wavevec- 
tors satisfy the excitation energy expression — = ujo for kp > 7r/4. For kp < 7r/4, only the region 
t/coo > l/[4sin(2A;F) sin(A;F)] is relevant as Rcxo(l5, wo) = lies above this line (see Fig. 3 of main text). For 
kp > 7r/4, let the two wavevectors that satisfy — eg = be fci and k2- We will use the blocking method to block 

the states ki—p and ^2— pby the two electrons on the Fermi surface. The state = \kF^ k\—p, —kp^ ^2— p; n-p) 
leads to the energy 

E^n = + Wo + 4i cos(A;f) + nwo[l + g^oJoRexoip, wo)] 

+g'^wlRe[Ilo{ki - p, kp,u!o) + Ilo{k2 - p, -kp, wq)]- (C8) 

Here also the macroscopic instability occurs when g> Qc with obtained from Eq. (C6). For extreme values of kp, 
i.e., kp close to 7r/4 or 7r/2, the large n state |$'') is the lowest energy state for g > g,; < gi (see Fig. 8). On the 
other hand, at intermediate values of kp, the state |$") is the stable state for gi < g < gc (see Fig. 8). 
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